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Abstract 
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our results for an investigation of the issue of zero mode degeneracy. 
We construct explicitly a class of multiple zero modes with their gauge 
potentials. 
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1 Introduction 



The Abelian Dirac equation in Euclidean space in two and three dimensions 
is 

= a k (-id k - G k (r))^(r) = 0, (1) 

where, in three dimensions, r = (x,y,zY, k = x,y,z, a k are the Pauli ma- 
trices, and G k = (G x , G y , G z ) is the gauge potential (we use superscripts for 
the components because subscripts will be exclusively reserved to describe 
partial derivatives). In two dimensions all z components are absent. Further, 
\1/ = ($, xY is a two-component spinor. 

The Pauli equation is obtained from the the Dirac equation by simply 
squaring the Dirac operator V acting on \1>, 

_ G k f - a l H l ) * = 0, (2) 

where H is the magnetic field, H = V x G. Solutions of the Dirac equation 
(1) are, at the same time, solutions of the Pauli equation, therefore we shall 
treat them on an equal footing. 

Solutions to the Pauli or Dirac equation are relevant in several instances. 
They describe the behaviour of non-relativistic electrons in the presence of 
magnetic fields, and their existence, among other issues, influences the sta- 
bility of non-relativistic matter [1], [2]. On the other hand, they influence 
the behaviour of the Fermion determinant det V for relativistic electrons and 
are, therefore, relevant for the strong field behaviour of relativistic electrons 
[3]-[5], and for a proper path-intergral quantization of QED in two [6]-[8] 
and three dimensions. 

For the Dirac equation in two dimensions the most important information 
on solutions is provided by the Aharanov-Casher theorem [9], which states 
that there are n square-integrable solutions when the magnetic flux divided 
by 2tt is between n and n + 1, i.e., n < (flux/27r) < n + 1. For the Dirac 
equation in three dimensions the available information is much scarcer. The 
first examples of solutions have been given in 1986, see [2]. Further examples 
have bee provided in [10], [11], [12], and the existence of multiple solutions 
(zero mode degeneracy) has been first demonstrated in [13], [14], and [15]. 
Further results on zero-mode supporting gauge potentials may be found in 
[16], [17], and in [18]. A general classification of zero-mode supporting gauge 
fields is still missing. 
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In this paper we study the symmetries of the Dirac equation (1) both 
in two and three dimensions, in order to gain some further insight into the 
solution space of these equations. In Section 2 we calculate the symmetries 
of the Dirac equation in two dimensions and briefly apply our results to the 
issue of multiple zero modes. In Section 3 we calculate the symmetries of 
the Dirac equation in three dimensions. We use the method of prolongations 
for all our symmetry calculations, which is described in detail in the book 
[19]. In Section 4, we use our results on symmetries of the three-dimensional 
Dirac equation for a discussion of the issue of zero mode degeneracy (i.e., 
multiple solutions) in three dimensions. In Subsection 4.1, we discuss how 
our results on symmetries can be used to construct multiple zero modes in 
principle. In Subsection 4.2, we apply the results of Subsection 4.1 for the 
explicit construction of a class of multiple zero modes together with their 
corresponding gauge potentials. 

2 Symmetries of the Dirac equation in two 
dimensions 

For the Dirac equation in two dimensions it is known from the Aharanov- 
Casher theorem that all zero modes (i.e., square-integrable solutions) are 
either left-handed (i.e., the lower component of ^ is zero) or right-handed 
(the upper component of \1/ is zero). Further, a solution of the first type 
(left-handed) may be mapped into a solution of the second type by the simple 
replacement G k — > — G k , therefore we may restrict, e.g., to the left-handed 
case 

(-id x - G x + d y - iG y )<$> = 0. (3) 

By introducing the modulus and phase of $ via In $ = p + iX we may rewrite 
Equation (3) in terms of two real first order equations as 

Ax = Py + X x -G x = (4) 
A 2 = Px -X y + G y = (5) 

where from now on subscripts denote partial derivatives, i.e., p x = d x p, etc. 

The vector field generating generic transformations on the independent 
and dependent variables is 

v = Xd x + Yd y + R8 P + Ld\ + F x d Gx + F y d Gy (6) 



2 



where X, Y, R, L, F x , F y may depend on all independent and dependent vari- 
ables. The equations (4, 5) contain first derivatives of the variables A and p, 
therefore we need the first prolongation of v w.r.t. these variables, 

pr v = v + R*d Px + R y d Py + L x d\ x + L y d Xy , (7) 

where, e.g. 

R x = D X R - p x D x X - p y D x Y. (8) 

Here D x is a total derivative, and we give the explicit (rather lengthy) ex- 
pressions in Appendix A. 

Now, a symmetry of a PDE F(p, pj, . . .) = (of n-th order, say) is a 
solution of the equation pr^v(F) = which holds on-shell, i.e., when the 
original PDE is used together with its prolongations (PDEs that follow from 
F = by applying total derivatives). As said, for the details of the used 
formalism we refer to the book ([19]). 

Concretely, the determination of the symmetry transformation proceeds 
as follows. We require that prt>(Ai) = = prt>(A 2 ) whenever the equations 
(4) and (5) hold. Explicitly this means that 

R y + L x -F x = (9) 
R x -L y + F y = (10) 

whenever the equations (4) and (5) hold. Equations (9) and (10) contain a 
number of algebraically independent functions (like p x , p x G y , etc.) multiply- 
ing the coefficients (like R x , R p , etc.) which we want to determine. We may 
use the equations (4) and (5) to eliminate, e.g., and \ y from the equations 
(9) and (10). Then we demand that the coefficient of each algebraically inde- 
pendent function (like p x , p x G y , etc.) vanishes separately. This leads to more 
than 30 conditions for each of the two equations (9) and (10). Fortunately, 
most of these conditions are quite trivial, like, e.g., X\ = = Xa*, etc. The 
final result is that Eq. (9) leads to the conditions 

X = X(x,y), Y = Y(x,y), (11) 

R = R(p,X,x,y) , L = L(p, A, x, y) , (12) 

R p -L x = Y y -X x (13) 
L p + R x = X y + Y x (14) 
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and 

F X = L X + R y + (L x - X X )G X + (R x - Y x )G y (15) 
whereas Eq. (10) leads to the conditions 

X = X(x,y), Y = Y(x,y), (16) 

R = R(p,X,x,y) , L = L(p, A, x, y) , (17) 

R P -L X = -Y y + X x (18) 
L p + R x = -X y -Y x (19) 

and 

F y = L y — R x — (R x + X y )G x + (L x - Y y )G y . (20) 
For later convenience we introduce the complex notation 

z = x + iy, Z = X + iY, a = p + iX , T, = R + iL (21) 

and 

G = ]-{G x + iG v ) , F = ]-{F x + iF v ) (22) 

which implies 

d c = \{d p -id x ) (23) 

etc. Compatibility between the two sets of equations requires 

X x = Y y X y = —Y x (24) 

which are just the Cauchy-Riemann equations, and 

L p = L x , R x = -L p , (25) 

which are the Cauchy-Riemann equations w.r.t. the target space variable 
a = p + iX. Therefore, X and Y are the real and imaginary part of a 
holomorphic function, 

X + iY = Z(z) , z = x + iy (26) 

whereas L and R are the real and imaginary part of a holomorphic function 
in the variable a, 

R + iL = Z(a,z,z), (27) 
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where E may still depend on z,z. Finally, for F we find 

F = ^(F* + iF") = -iE 2 - - GZ 2 - + GE CT . (28) 

The most general symmetry generator 

u = Z<9 2 + Z<9 2 - + E<9 CT + Ec^ + F8 G + F0 e (29) 

may be expressed as the semi-direct sum of a symmetry generator vz w.r.t. 
the holomorphic function Z(z) and a generator t> s w.r.t. the function E(<7, z, z), 
where 

v z = Zd z + Z9 2 - GZ 2 -«9 G - GZ z dc (30) 

and 

t; s = Ed a + td 9 + HE 2 + GE CT )«9 G + (*E 2 + G%)d G . (31) 
Further, they obey the infinite-dimensional Lie algebra 

[vznVzz] = v Za , Z 3 = Z X Z' 2 - Z 2 Z[ , (32) 

[^Si , vy, 2 ] = v-s 3 , E 3 = E 1 E 2)(7 - E 2 E 1;(7 , (33) 

which are just two copies of the Virasoro algebra in coordinate and target 
space, and 

[v z , v^]=v t , E = E 2 Z + E 2 -Z . (34) 

It is quite interesting that the symmetry turns out to be so large, the semi- 
direct product of two conformal groups in two dimensions. 

[Remark: the symmetry transformations we found cover the whole solu- 
tion space in the sense that any solution $ of the Dirac equation (3) for any 
gauge potential G can be found by applying a symmetry transformation to 
the trivial solution $ = 1, G — 0. In fact, a target space transformation E 
of the type E = f(z,z) is sufficient. For this sub-class of transformations 
the above Lie algebra is Abelian, therefore the exponentiation for finite E is 
trivial.] 

As an application let us briefly discuss the issue of degeneracy of zero 
modes (i.e., multiple solutions of the Dirac equation for one gauge field G). 
The condition that the gauge field does not change is, of course, that the 
F in Eq. (28) is zero, which leads to the conditions S CT = S 2 = 0, that is 
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E = g(z) is a function in the variable z only. As said, the exponentiation of 
such £ is trivial, therefore a 

$' = e 9(z) $ (35) 

is a local zero mode for the same gauge potential as $. Single- valuedness of 
$' on the whole Euclidean plane restricts the allowed g (expg(z) should be 
single-valued) and the condition of square-integrability further restricts the 
allowed g(z). The fact that further zero modes of the same Dirac operator 
may be produced by multiplying the given one with analytic functions of z 
is, of course, well-known. 



3 Symmetries of the Dirac equation in three 
dimensions 

Introducing the real and imaginary part of the spinor components via $ = 
a + i/3, x — 7 + 2<5, the Dirac equation (1) in three dimensions is equivalent 
to the four real equations 

/ = f3 z + 5 X - ly - G z a - G x 1 - G y 5 = 

// = a z + 7 x + 5 y + G Z P + G x 5 - G Vr ) = 

/// ee P x + a y — 8 Z — G x a + G y /3 + G z ^ = 

IV ee a x - P v - 7 2 + G x /3 + G y a - G z 5 = 0. (36) 

If viewed as a system of algebraic equations for the three real unknowns 
(G x ,G y ,G z ), these four equations must be algebraically dependent. The 
dependence is given by the condition 

V ee V • E = 2(a 7 + pS) x + 2(a5 - p>y) y + (a 2 + (3 2 - 7 2 - 5 2 ) z = (37) 

where 

E ee 



( 2(af + f35) 

2{aS - Pi) \ (38) 



V a 2 + [3 2 - 1 2 -5 2 



is the spin density of the spinor 
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Again, we want to study the symmetries of the above equations (36) using 
the method of prolongations. The symmetry generating vector field is 

v = xd x +Yd y +zd z +Ad a +Bdp+cd 1 +Dd s +F x d Gx +F y d G y+F z d G , (39) 

and its prolongation to the needed order is 

W v = v + A x d ax +A y d ay +A z d az +B x d f3x + B y d (3y + B z d Pz 

+C x d lx + C y d ly + C z d lz + D x d Sx + D y d 5y + D z d 5z (40) 

where, e.g., A x is given by 

A x = A x + A a a x + Ap(3 x + A^ x + A S 8 X - a x X x - a y Y x - a z Z x (41) 

Here, we have already used some first results of the calculation below, namely 
that (analogously to the case in two dimensions) no coefficient depends on 
the gauge potential, and that the coefficients X, Y, Z of the translations in 
base space only depend on the base space coordinates x, y, z (otherwise the 
resulting expression (41) would be more complicated, analogously to the R x 
of the two-dimensional case, which is given in Appendix A). 

Acting with the prolonged vector field on the four equations (36) (pr v(I) = 
0, etc.) leads to the following set of four equations 

G X C + G y D + G Z A + F x 1 + F v 5 + F z a - D x + C y - B z = 
G x D-G y C + G z B + F x 5-F y 1 + F z (3 + C x + D y + A z = 
-G x A + G y B + G z C-F x a + F y l3 + F z ~i + B x + A y -D z = 
G X B + G y A — G Z D + F x /3 + F y a — F z 5 + A x — B y — C z = 

(42) 

Next we have to insert the explicit expressions for A x etc. We obtain algebraic 
expressions in terms of the target space variables and their partial derivatives, 
where we may eliminate four derivative terms with the help of the Dirac 
equation (36). (Explicitly, we eliminated 5 x ,8 y ,S z and 7 2 .) For the resulting 
expressions one has to require that the coefficient multiplying each partial 
derivative of the target space variables (like a x or G' x ) or product of partial 
derivatives vanishes separately. Doing this one realizes quickly that nothing 
may depend on the gauge potential and that the coefficients X, Y, Z may 
only depend on x, y, z. 
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The remaining coefficients which are multiplied by at least one partial 
derivative of a target space variable serve to determine the coefficients X, Y, Z 
and A, B, C, D. They are of two types. The first type consists of pairs of 
equations like 

As + D a = Y z + Z y , A s + D a = -Y z -Z y ^A s + D a = 0, Y z + Z y = 

which lead to equations for the X, Y, Z and A, B, C, D independently. The 
equations for X, Y, Z are 

-X-X Yy Z Z 

X y + Y x = 0, X z + Z x = 0, Y z + Z y = (43) 

and have the conformal transformations in three-dimensional space as general 
solutions, 

Ql 

X = — (x 2 -y 2 -z 2 ) + 6 2 xy + 6 3 xz-l :i y + l 2 z + ax + x 
Y = — (y 2 -x 2 -z 2 ) + 6 1 xy + 6 3 yz + l 3 x-l 1 z + ay + y° 

n3 

Z = —(z 2 -x 2 -y 2 ) + e l xz + e 2 yz-l 2 x + l l y + az + z° (44) 

Here 9 l ,l\a and r° are constant parameters which parametrize the infinites- 
imal conformal transformations. 
The equations for A,B,C,D are 

A a = Bp = C 7 = Ds 

Ap + B a = 0, A 7 + C a = , A s + D a = 

B^ + C p = 0, B 5 + Dp = 0, C s + Dt = Q (45) 

and have the conformal transformations in the four- dimensional target space 
as solutions, 

A = ^(a 2 - (3 2 - ^ 2 - 5 2 ) + ( 2 af3 + ( 3 ai + ( 4 a5 + 

Xa - v 3 [3 + w 2 7 - v l 5 + a 
B = ^([3 2 -a 2 - 1 2 -5 2 ) + ( 1 a[3 + ( 3 [3 1 + ( 4 [36 + 



8 



\(3 + v 3 a - - u 2 5 + (3° 

C = ^( 7 2 -« 2 -/? 2 -5 2 ) + C 1 «7 + C 2 /37 + CV + 
A7 - v 2 a + v}(3 - u 3 5 + 7 

D = ^(5 2 -a 2 = /3 2 -7 2 ) + CW + C 2 ^ + CV + 

X5 + v l a + u 2 [3 + m 3 7 + 5° (46) 

Here the parameters ( j ,u\v l ,\ and a , (3°, 7 , 5° may still depend on the 
base space coordinates (x,y,z). 

The second type of equations consists of equations like C$ — Ap = X y and 
establishes relations between the A, B, C, D and the X, Y, Z, thereby further 
restricting the possible A, B, C, D. The result is that the proper conformal 
transformation on target space must be absent, C J = 0, and that the six 
target space rotation parameters are no longer independent, 

^-v 1 = I 1 + 2 z — 6 3 y = L 1 

U *- V * = l 2 + p x _ d l z = L 2 

u 3 -v 3 = I 3 + 6 1 y — 6 2 x = L 3 (47) 

leading to 

A = \a-v 3 [3 + w 2 7 - v 1 5 + a 

B = X(3 + v 3 a - {v 1 + L x ) 7 - {v 2 + L 2 )5 + (3° 

C = \ 1 -v 2 a + (v 1 + L 1 )(3-{v 3 + L 3 )5 + 1 ° 

D = X5 + v 1 a + {v 2 + L 2 )(3 + {v 3 + L 3 ) 7 + 5° (48) 

where A, v l and a , . . . may still depend on (x,y,z). By shifting v 1 — > v l — 
(L l /2) the L l are distributed more symmetrically, 



A 


= Aa - (w 3 - 


r 3 


(, 2 - 


^)7- 
2 ; ' 


(v 1 - 


Y )6 + a ° 


B 


= A/5 + (w 3 - 


Y )a ~ 


(^ + 




(v 2 + 




C 


= A 7 -(t; 2 - 


Y )a + 


(^ + 




(v 3 + 


|W 


D 


= A5 + (v 1 - 


y)« + 


(t; 2 + 


T 2 


(v 3 + 


^)7 + ^ 



(49) 
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and we find the half-integer valued representation of the base space rota- 
tions (with parameters l z ) acting on the (spinor) target space variables. It 
is interesting to note that the half-integer values appear also for the pa- 
rameters of the proper conformal transformations on base space (remember 
U = I' 1 + e ljk 9 j r k ). 

[Remark: for a Dirac equation without the additional condition V- £ = 0, 
see Eq. (37), the above expressions would be the final result. This is the 
case, for instance, for the Dirac equation for a Weyl (i.e., two-component) 
spinor in Euclidean space in four dimensions, for which expressions analo- 
gous to (49) would constitute the final result (with the conformal base space 
transformations in 4 instead of 3 dimensions, of course).] 

We have used all information from the coefficients of equations (42) which 
contain at least one partial derivative of a target space variable. It remains 
to evaluate the parts without a partial derivative (i.e., the coefficients of the 
identity). These contain the coefficients F x , F y , F z linearly and are therefore 
just the determining linear equations for these coefficients. However, they 
give 4 equations for 3 unkonwns, therefore they must be linearly dependent, 
which leads to one further constraint equation. An easier way to calculate 
this constraint is to calculate the action of the prolonged vector field (40) on 
the original constraint (37), 

prv(V) = (50) 

or explicitly 

A x -f + a x C + A lx + aC x + B x 5 + (3 X D + B5 X + (3D X + 

A y 5 + a y D + A5 y + aD y - B y 1 - (3 y C - B ly - (3C y + 

A z a + a z A + B z {3 + (3 Z B - C 2 7 - lz C - D z 5 -5 Z D = (51) 

Inserting the explicit expressions for A x , etc. leads to an expression con- 
taining both coefficients multiplying partial derivatives of the target space 
variables (like a x ) and a coefficient of the identity 1. After eliminating four 
partial derivatives (e.g. S x ,5 y ,5 z and j z ) with the help of the e.o.m., one 
realizes that the coefficients of the remaining partial derivatives vanish iden- 
tically. It remains to evaluate the coefficient of the identity. This leads in 
fact to four conditions, because the use of the e.o.m. has reintroduced the 
gauge potential functions G x ,G y , G z into the above expression. As none of 
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the coefficients (A, etc.) may depend on the gauge potential, the total co- 
efficient multiplying each gauge potential function has to vanish separately, 
as well as the remainder. The conditions that the coefficients mulitplying 
G x ,G y ,G z vanish leads to 

a° = /3 = 7 = 5° = 0, v l =v 2 = Q (52) 

and the condition that the remainder vanishes leads to 

(VA + #)-£ = (53) 

With the shift 

a^a-0, e = e-f=e 1 x + e 2 y + e 3 z (54) 

and the definition 

= -v 3 (55) 

we therefore arrive at 



A = 


\a + (fi/3 


r3 

-Qa + —[3 
2 


L 2 


2 




B = 


A/3 — 4>a 


L 3 

-Q(3 a 

1 2 


L 1 

2 7 


L -6 
2 




C = 


A7 + <p5 - 


L 2 

- ©7 + Y a ' 




*s 

2 




D = 


\5 — 07 


-Q5-—a- 
2 


r 2 

4- 9+ 


T 7 


(56) 



with 

(VA) • E = 0. (57) 

Here A and are scale and gauge transformations on (spinor) target space, 
whereas the remainder is the representation on spinor space of the infinites- 
imal base space symmetries. 

However, there is a problem with the constraint (57). The constrained 

— # 

function A is still algebraically independent of £ (i.e., of the a, (3, etc.), 
therefore the constraint does not contradict our basic assumptions. The 
problem is that the Lie algebra of all infinitesimal transformations does not 
close on the constraint (57). E.g., the commutator of a target space scale 
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transformation and a base space rotation about the z-axis leads to another 
target space scale transformation which no longer obeys the constraint (57) 

[v x ,v P ]=v- x , (VA)-E^O (58) 

This is, in fact, not a surprise, because the gradient of the rotated scale 

function A must be perpendicular to the rotated vector E, (VA) ■ X = 0, and 
analogously for other base space transformations. 

Therefore, in order to have a closing Lie algebra of base space and target 
space transformations, we have to restrict to constant target space scale 
transformations, A = const. 

However, as transformations with non-constant A obeying the constraint 
(57) map solutions of the Dirac equation to new solutions and are of some 
interest as such, one may instead choose another restriction by restricting to 
the target space transformations parametrized by A and (i.e., by setting 
all base space transformation parameters equal to zero). The resulting Lie 
algebra is abelian and closes therefore trivially. 

We still have to calculate the coefficients F x , F y , F z from the coefficients 
of the identity of Eqs. (42). After a lenghty calculation one obtains the 
simple result 

F = -(Q + a)G + L x G- V0 - 4-S x VA (59) 

I ^ I 

4 Multiple zero modes 

In this section, we would like to apply the results on symmetries of Section 3 
to the issue of zero mode degeneracy (i.e., multiple solutions for a Dirac equa- 
tion with the same gauge potential). We shall give a more general discussion 
in Subsection 4.1, whereas we will provide some explicit, new examples of 
multiple zero modes in Subsection 4.2. 

4.1 General discussion 

The condition that a target space 
not change the gauge field is F = 

V0 = 



symmetry transformation of a spinor does 
0, see Eq. (59), or 

= — l -t x VA (60) 
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which, together with condition (57), implies that S, V0 and VA are mutually 
perpendicular. Further, Darboux's theorem tells us that (at least locally) we 

— * 

may express S like 

s = ve x ve (6i) 

where £ a = £ a (r), a = 1,2. If we now assume that and A are functions 

— * 

of £ a only, then their gradients are automatically perpendicular to E, and it 
remains to solve 

|V0| = |VA| , V0-VA = O. (62) 

This problem can be solved relatively easily for a subclass of functions £ a 
which fulfill one additional requirement, namely that the scalar products of 
the gradients of the £ a can be expressed in terms of the £ a again, up to a 
common factor, that is 

Vr-V£ 6 = Mr)<7 ab (0- (63) 

Here g ab plays the role of a metric in some two-dimensional space parametrized 
by the coordinates £ a . If Eq. (63) holds, then Eq. (60) simplifies to 

4> a e ac = g cb \ b (64) 

where a = d^ a (f), etc. and 

~ e ab = g \ e "\ g = det(g ab ) = g n g 2 2 - #i 2 #2i (65) 

and e ab is the usual antisymmetric symbol in 2 dimensions. Eqs. (64) can be 
solved by observing that they just provide a generalization of the Cauchy- 
Riemann equations to the case of a general surface. We only have to find the 
coordinate transformation from £ a to some new coordinates £ /a such that the 
metric g' ah w.r.t. the new coordinates is conformally equivalent to the flat 
metric 5 ab , that is 

3' ab = ^ d £l9 cd = KOt cd - (66) 



This problem always has a solution in two dimensions [21]. In terms of the 
coordinates £ /a Eqs. (64) read 



<p a e ac = 5 cb \ b (67) 
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which are just the Cauchy-Riemann equations. Therefore, our problem is 
solved by an arbitrary complex function u of the complex variable (' where 

«(C0=exp(A + #), C^^ + tf. (68) 

Here "solution" means that if a spinor \I/ with = E solves a certain 

Dirac equation, then u^/ locally solves the same Dirac equation, 

= => Vv& = 0. (69) 

Regularity requirements further restrict the allowed u. For instance, the 
functions £ a need not be well-defined in all IR 3 (they usually are not), but 
u certainly has to be well-defined (which may not be possible, in which case 
the problem has no acceptable solution). The condition that is square- 
integrable restricts to a finite number (which, again, may be zero) of linearly 
independent functions u. It should be emphasized that all known examples 
of multiple zero modes are of the above type (69), see [13]— [15]. It is an 
interesting open question whether there exist other types, as well. 

Finally, let us point out that the procedure described above can certainly 
be reversed. That is to say, choose a pair of functions £ a which obey Eq. (63) 
and calculate the corresponding £ /a via Eq. (66). Then u^> with u given by 
Eq. (68) will provide additional local zero modes for a whole class of spinors 
\1/ with spin densities given by E = F(^ a )V^ 1 x V£ 2 for (almost) arbitrary 
real functions f(i a ) (of course, / should be well-defined in all IR 3 ). It is not 

— * 

difficult to reconstruct the spinor \1/ from the spin density E and the gauge 
potential G from the spinor \P (see, e.g., [2]). We shall explicitly demonstrate 
how this works by constructing a class of multiple zero modes in the next 
subsection. 

[Remark: Equations (62) can be expressed in terms of the complex vari- 
able u = exp(A + i(f>) like 

(Vw) 2 = (70) 

which is the complex static eikonal equation. Further, condition (63) is 
equivalent to the condition that u - when interpreted as a map from one-point 
compactified IRq to one-point compactified (D - is a Riemannian submersion 
up to Weyl transformations (local conformal rescalings of the metric). That 
is to say, u is a composition of maps 

u '. IR 3 ^ M 3 R 4- M 2 ^ (Co (71) 



14 



where W is a Weyl transformation, R.S. is a Riemannian submersion, and 
Ai 3 and M 2 are compact manifolds in three and two dimensions, respectively. 
A detailed discussion of these issues can be found in [20]. In [15], Rieman- 
nian submersions were used to construct multiple zero modes within a more 
geometrical context.] 



4.2 A class of multiple zero modes 

Here we want to construct explicitly a class of multiple zero modes by starting 
with a pair of functions £ a which obey Eq. (63), as explained in the last 
subsection. In particular, we want to make use of the results of [20], where a 
class of Hopf maps obeying the eikonal equation (70) were found, and where 
the geometric explanation for their existence was provided. In concordance 
with these results, we therefore choose the functions 

= In sinh 77 = - InT, £ 2 = md + nip (72) 

where we introduced toroidal coordinates (77, ip) (and, for later convenience, 
the variable T = sinh 2 77) related to the cartesian coordinates (x, y, z) via 

t- • 1,2 4Qr 2 + y 2 ) 2z y 

1 = sinh 77 = — r — , v = arctan — , ip = arctan — . (73) 

Az 2 + (r 2 - l) 2 r 2 - 1 x 

Further, m and n are nonzero integers such that exp(27ri/£ 2 ) is a single- valued 
function for integer I (the geometric significance of the integers m and n is 
that they provide the complex function (Hopf map) / = exp^ 1 + i^ 2 ) with 
the Hopf index H = mn). 

Following the results of [20], it can be shown that the pair £ x ,£ 2 obeys 
Eq. (63), and that a new pair £ /2 obeying Eq. (66) can be found without 
difficulty Indeed, one finds easily that £ a obey Eq. (63) with 

, . . (cosh 77 — cost?) 2 cosh 2 77 . , 

h(r) = — 2 (74) 

sinh 77 

11 1 12 21 n 22 m 2 sinh 2 r] + n 2 

9 =1, 9 =9 =0, 9 = T2 • ( 75 ) 

cosh 77 

Hence, the induced metric g ab is already diagonal but not yet conformally flat. 
However, as g 22 only depends on (i.e., on rf), a coordinate transformation 
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involving only £ x is sufficient, £ x — > £ /2 = £ 2 , such that g 11 — > g' 11 = 

g 22 = g /22 . Explicitly, the transformation reads 



{" = In 



/ / / \ M \ 

( Iml coshr? + \/n 2 + m 2 sinh 2 "■ ! 
sinhH ' ^ ' ' 



rj 



V 



n| cosh 77 + \Jn 2 + m 2 sinh 2 77 



|n| 



(76) 



see [20]. The function (' = C,' 1 + i£ /2 itself is not single- valued, but the 
function 

C = expC' = exp(e ,1 + ^ 2 ) (77) 

is. It follows that for each spinor \I/ with spin density E = ty^aty = F(C, a ) V^ 1 x 
V£ 2 (which is a formal zero mode for some gauge potential, because S obeys 

— * 

V-S = 0), u(()^/ are further formal zero modes for the same gauge potential, 
where u is a rational function of its argument. A customary basis for the 
functions u is u — for integer I. 

The formal zero modes described so far (implicitly via their spin den- 
sity) are single- valued functions on all IR 3 , but we have not yet taken into 
account the condition of square-integrability. Before doing so, we want to 
make some simplifying assumptions on the type of spin density we want to 
discuss. Firstly, we assume that the function F depends on (i.e., on 77 or 
T) only. Secondly, we reexpress F like F = e M , that is, we write for a 
general spin density 

g(M) = e M(T)g(o) ^ g ( o) = * T Ve 1 x V£ 2 . (78) 

The reason for this is that is a well-behaved and integrable spin density, 
whereas V^ 1 x V£ 2 is not (clearly, integrability of the spin density / rf 3 r|S| < 
00 is the same as square-integrability of the corresponding spinor). Another 
reason for the choice of E^, which does not concern us much here, is the fact 
that S(°) is the Hopf curvature for the Hopf map exp^ 1 + i^ 2 ). In cartesian 
coordinates, reads 

j~ / 2nxz — 2my 
= - | 2 ^ 3 I 2nyz + 2mx 1 (79) 
"~ V ' ! n(l-r 2 + 2^ 2 ) 



:i+r 2 ) 



16 



and its spinor *(°) with = lr(o)t^r(o) is 

*<°> = 2 ^ (80) 

(1 + r 2 )f ^(1 + r 2 )r + n(2z 2 + 1 - r 2 ) 



where 



' n(2^ 2 + 1 - r 2 ) + (1 + r 2 )r \ 
2(x + iy)(nz + im) J 

n 2 + m 2 T \ 5 f n 2 + m 2 sinh 2 77 " 



(81) 



(82) 



\ 1 + T J \ cosh 2 77 / 

and we chose the gauge such that the upper component is real. The calcu- 
lation of the spinor from the spin density is explained in Appendix 
B. The spinor (81) is regular everywhere. Further, it is a zero mode for some 
gauge potential G^°\ by construction (this gauge potential is well-behaving 
and leads to a square-integrable mag netic field i? (0) = V x G (0) ; its explicit 
expression, which is quite lengthy, is displayed in Appendix B, together with 
an explanation of its calculation). Additional formal zero modes for 
the same gauge potential are not square-integrable, i.e., the corresponding 
spin density |£*£|*E(°) is not integrable (for either positive or negative integer 
I), as may be inferred easily from the explicit expression for and from 
the limiting behaviour 

Km|C*CI~rW, lim|C*CI~rM. (83) 

1 1 -^00 

But with an appropriate choice of M, it is easy to find spin densities E( M ) 
such that |C*C|'E( M ) remains integrable for some non-zero values of I. 

Concretely, we assume that M(T) and M'(T) are finite for all finite values 
of T, i.e., 

\M(T)\< 00 A \M'(T)\< 00 for T < 00, (84) 
and that M behaves for large T like 

lim M(T) mk In T + M(T) (85) 

where the remainder M has to obey 

lim \TM(T)\ < 00. (86) 

T—>oo 
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Here the conditions (84) and (86) are chosen in order to avoid both spurious 
(i.e., pure gauge) and physical singularities for the corresponding gauge po- 
tential, and condition (85) is chosen in order to have a non-trivial result (i.e., 
multiple square-integrable zero modes). Condition (85) induces a spurious 
(gauge) singularity in the corresponding gauge potential at T = oo which 
has to be cured by an appropriate gauge fixing. Consequently, the spinor 
\$( M ) f or the spin density X/ M ) with the appropriate gauge fixing is 

where the additional factor exp(ikm , &) provides the gauge fixing, see below. 
Under the above assumptions the spinors 

^(M) = C ^(M) ? l = o... k (88) 

are all square-integrable zero modes for the same gauge potential & M \ 
It remains to determine G^ explicitly, which will be obtained with Eq. 

(59) of Section 3. In fact, using Eq. (59) (where M/2 ~ A and kmd ~ 0; 
remember that Eq. (59) also holds for finite A and 4> because of the abelian 
nature of the target space symmetry transformations) one easily calculates 

M'T 

q(m) = q(o) _ __( mW + nVv? ) _ mA;W ( 89 ) 

(as said, the expression for is provided in Appendix B). Here the addi- 
tional pure gauge term precisely cancels a pure gauge singularity at T = oo, 
as may be checked easily. Further, we may see how it works that multiplica- 
tion of the spinor \1/( M ) by £ does not change the gauge potential. The crucial 
point is that multiplication by ( corresponds to choosing M = |C*CI> an d for 
this M it holds that ^jL = —1, so that the M dependent term in (89) is, 
in fact, pure gauge and is cancelled by the contribution Varg£ = V£ 2 , see 
(72). This just shows explicitly that M/2 = f 1 ~ A and £ 2 ~ <p fulfill Eq. 

(60) which, of course, must be true by construction. 

Therefore, we have succeeded in constructing explicitly a class of multiple 
zero modes together with their gauge potentials starting from the functions 
(72) (or, equivalently, starting from the higher toroidal Hopf maps of Ref. 
[20]). For the simplest Hopf map (i.e., for m — n — 1) these zero modes have 
already been obtained in [13, 14, 15], whereas for the higher toroidal Hopf 
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maps they are new. 



Appendix A 

The coefficients for the prolongation prt> of the vector field v in Section 
2 are 

R x = D x R-p x D x X- Py D x Y (90) 

W = D y R — p x D y X — PyDyY (91) 

L x = D X L — \ X D X X — \ y D x Y (92) 

L y = D y L — \ x DyX — XyDyY (93) 

where D x and D y are total derivatives w.r.t. x and y, respectively. Explicitly 
R x reads 

R x = Rx + RpPx + R\^x + Rg*G x + RgvG x 

- Px (X x + X pPx + X x \ x + X G *G x x + X G yGl) 

Py (Y x + Y pPx + Y x \ x + Y G *G x + Y G yGl) (94) 

and analogously for the other coefficients. 



Appendix B 

The calculation of the spinor from the spin density proceeds as 
follows. For a gauge such that the upper component of the spinor is real, a 

— * 

spinor \l> may be obtained from its spin density S formally by the algebraic 
relation, see Ref. [2], 

* = I - { Tit ? 5 ) • (95) 

— * 

This gauge is globally admissible provided that the third component of S is 
positive whenever its first and second component are zero, i.e., 

S 1 = A S 2 = =>• S 3 > (96) 
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which holds for the spin density Y*^ of Eq. (79). Therefore, the spinor 
of Eq. (80) may be computed with the help of the above expression. 

The calculation of the gauge potential for the spinor (zero mode) 
proceeds as follows. Remember that once a zero mode * is given, the Dirac 
equation may be viewed as an overconstrained system of linear, algebraic 
equations for the three components of the gauge potential. If the constraint 
is fulfilled, this system can be solved explicitly and results in the following 
expression for the corresponding gauge potential G, see Ref. [2], 

G= -U (V x S + 2Im(^ t W)) . (97) 
2 j XI I 

Inserting the explicit expression for into the above formula leads to the 
following rather lengthy expression for G^°\ 



—n(5 — r 2 )y + 6mxz ^ 
^ i+r > \ m(Qz 2 + 2-4r 2 ) J 



+ (98) 



T(l + r 2 )[r(l + r 2 ) + n(2z 2 + 1 - r 2 )] 



/ — (n 2 z 2 + m 2 )y \ 

(n 2 z 2 + m 2 )x 
\ mn(r 2 - z 2 ) ) 



(99) 



(for a direct comparison with the results of [13, 14] for the case m — n — 1, 
one must take into account the fact that the gauge chosen in [13, 14] differs 
from the one chosen here by the gauge function <fi = arctan(,z)). Further, 
Eq. (97) shows that G is singular whenever |E| is zero. This singularity 
may be a spurious (gauge) singularity, in which case it can be cured by an 
appropriate gauge fixing as in Eq. (89), or it may be physical, in which case 
the corresponding zero mode is not admissible (of course, all the zero modes 
given in Subsection 4.2 are admissible). 
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